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ABSTRACT 

We investigate the (large-scale) bar fraction in a mass-complete sample of M > 1O 1O ' 5 M0 
disk galaxies at 0.2 < z < 0.6 in the COSMOS field. The fraction of barred disks strongly 
depends on mass, disk morphology, and specific star formation rate (SSFR). At intermediate 
stellar mass (10 10 5 < M < 10 11 M@) the bar fraction in early-type disks is much higher, 
at all redshifts, by a factor ~2, than that in late-type disks. This trend is reversed at higher 
stellar mass (M > 10 11 M©), where the fraction of bars in early-type disks becomes sig- 
nificantly lower, at all redshifts, than that in late-type disks. The bar fractions for galaxies 
with low and high SSFRs closely follow those of the morphologically-selected early-type and 
late-type populations, respectively. This indicates a close correspondence between morphol- 
ogy and SSFR in disk galaxies at these earlier epochs. Interestingly, the total bar fraction in 
10 10 5 < M < 10 11 Mq disks is built up by a factor of ^2 over the redshift interval explored, 
while for M > 10 11 M Q disks it remains roughly constant. This indicates that, already by 
z ~ 0.6, spectral and morphological transformations in the most massive disk galaxies have 
largely converged to the familiar Hubble sequence that we observe in the local Universe, while 
for intermediate mass disks this convergence is ongoing until at least z ~ 0.2. Moreover, these 
results highlight the importance of employing mass-limited samples for quantifying the evo- 
lution of barred galaxies. Finally, the evolution of the barred galaxy populations investigated 
does not depend on the large-scale environmental density (at least, on the scales which can be 
probed with the available photometric redshifts). 

Key words: Galaxies - galaxies: fundamental parameters - galaxies: bars - galaxies: forma- 
tion. 



1 INTRODUCTION 



Numerous observational studies have demonstrated that large- 
scale stellar b ars are remarkably common amongst local disk 
galaxies (e.g. de Vaucouleurs et al, 1991 1_ Eskridge et alj l200d : 
iMenendez-Delmestre et al JI2007I : Aguerri etlril 20091) . and that a 



substantial population of barred d i sks exists out to at least redshift 
unity (e.g. I Abr aham et al ] 1 19991 : Elmegreen, Elmegreen, & Hirst 
l2004IJogee et alj|2004 lElmegreen et alj|2007l ; ISheth et altoolT 
Quantifying the fraction of barred disks as a function of 
redshift, and its dependence on fundamental galaxy proper- 
ties, is an essential step towards understanding galaxy for- 



© 0000 RAS 



2 Cameron et al. 



mation in a cosmological context. Importantly, large-scale 
bars serve as signposts of m assive, dynamically-cold disks 
dAthanassoula & Seilwoodl ll986). constraining their epoch of for- 
mation dShefh et al.ll2008f) . Moreover, bars are key drivers of sec- 
ular evolution in their host galaxies, redistributing angular mo- 
mentum, enha ncing nuclear star formation , and building pseu 



2 DATASET AND ANCILLARY MEASUREMENTS 



dobulges (e.g. [Linden-Bell & Kalnaisl 19721; Combes et alj 



1990 



2002 



Knapen et all 19951; ICarolloll 19991; ICarollo et al.ll 19911200 ll. 
2007 1; ISheth et alJl200d ; lKormendv & KenmcuttbopllSheth et alJ 
2005; Deb attista et alj 12006 ; iFovle et alj 120091 ; IComeron et al.1 
20icl). 



A number of recent studies have examined the disk galaxy 
bar fraction for luminosity-s elected samples at intermediate-to - 
high redshif ts in the field (e.g.ljogee et all2004l;ISheth et al.l2008l) , 
in clusters dMarinova et aljfeooj^ . or both dBarazza et alj|2009l) . 
I Jogee et all d2004l) recovered a constant bar fraction out to red- 
shift unity in the GEMS survey for disk galaxies at Mv ^ —19.3 
and —20.6 mag. Conversely, for galaxies at L v (M v = —21.7 
mag at z = 0.9) and brighter in the Cosmological Evolution 
Survey (COSMOS), ISheth et al.1 d2008h i dentified a build up of 
the barred population ove r the same epoch. iMarinova et al. (2009) 
and iBarazza et al. I d2009h have identified a minimal dependence 
of the bar fraction on environment (see also iMendez-Abreu et al.1 
2010), except perhaps in cluster centres, where bar formation ap- 
pears to be enhanced. Both studies also recovered a rise in the 
bar fraction towards later morphologi cal types, or blu er, increas- 
ingly disk-dominated systems, while ISheth et al.1 d20080 recovered 
an enhanced bar fraction in redder, increasingly bulge-dominated 
systems at high redshifts. However, care must be taken when com- 
paring Jhejsej^sultsdueto the different luminosity limits adopted 
(e.g., IMarinova et alj|2009l and IBarazza et alj|2009l sampled faint 
disks down to M v ^ -18 mag at z ~ 0.165, and M v < -20 
mag at 0.4 < z < 0. 8, respectively) . Inter estingly, the high- 
est bar fraction found by IMarinova et al.1 ( 120090 was 75±11% for 
clumpy disks with distinct bulges at high luminosities (—21 < 
My < —20 mag). A further complication is introduced by the 
various methods of dis k selection employ ed, whether SED-fit type 
and visual inspectio n dShethetalJl2008f) , visual inspection alone 
dBarazza et a l. 2009), or colour , global Sersic index, and visual in- 
spection dMarinova et alj|2009l) . 

In the light of recent evidence that galaxy evolutionary histo- 
ries are tied closely to to t al stellar mass fe undv. Ellis, & C onselicel 
20051; iBaldry et alJl2006l ; lBolzonella et alJl2009l ; lTasca et alj|2009l ; 
Kovac et alj|2010h . the further study of evolution in mass-limited 



samples is a crucial step forward to undertake. In this Letter we 
quantify the bar fraction at intermediate-to-high redshifts, at fixed 
stellar mass, using a sample of 916 morphologically-classified disk 
galaxies in the COSMOS field. Specifically, we study the depen- 
dence of the (large-scale) bar fraction on detailed morphological 
type (early- or late-type disks), specific star formation rate (SSFR), 
and large-scale environmental density. 

The outline of this Letter is as follows. In Section [2] we re- 
view the COSMOS dataset employed in this work, and describe our 
sample selection. In Section|3]we explain our bar detection proce- 
dure and method of accounting for selection biases. In Section|4]we 
present our results, and in Section [5] we discuss their implications 
for galaxy formation scenarios. The construction and bar fraction 
analysis of a complementary sample of local Universe disk galaxies 
in the SDSS is described in the Appendix. We adopt a cosmolog- 
ical model with S1a = 0.75, Qm = 0.25, and h = 0.7, and all 
magnitudes are quoted in the AB system throughout. 



We use the COSMOS dataset dScoville et~al]|2007l) , consisting of 
ground-based and space-based, multi-wavelength imaging of an 
~2 deg 2 equatorial field. In particular, as the starting point for 
our sample definition, we adopt the /-band source catalogue of 
iLeauthaud etafl d2007l) extracted from th e COSMOS ACS /-band 
imaging frames dKoekemoer et alj 120071) , consisting of 156,748 
sources limited down to /ab = 26.6 mag. Careful artifact and 
star removal was performed to improve the robustness of this cata- 
logue, and computation of photometric redshifts was attempted for 
111,141 of these sources flagged as galaxies (i.e., removing stars 
and junk). 

Eleven photometric bands covering the wavelength range be- 
tween the COSMOS CFHT u and the Spitzer 4.5 band (see 
ICapak et al.ll2007l ; ISanders et alj|2007h were used to derive these 
photometric redshifts using our own Zurich Extragalactic B ayesian 
Redshift Analyzer (ZEBRaQ) code dFeldmann etai]|2006h . To this 
end, we used galaxi es with secure spectros copic redshifts from the 
zCOSMOS survey dLillv et al.ll2007l.l2009l) and modified with ZE- 
BRA our set of empirical templates dColeman. Wu. & Weedmanl 
1980; iKinnev et alj Il99 6) to account for systematic, template- 
dependent mismatches with the photometric spectral energy dis- 
tributions (SEDs) of the sample galaxies. An acceptable template 
match was able to be identified for all but 1695 of the 111,141 
input sources (a 1.5% failure rate). The resulting photometric red- 
shifts have an uncertainty of A(z)/(1 + z) ~ 0.023(1 + z) down 
to / < 22.5 mag, as directly tested with available spectroscopic 
redshifts from zCOSMOS. We estimated the photometric redshift 
accuracy at magnitudes fainter than / = 22.5 mag by dimming 
the photometry of zCOSMOS galaxies down to / = 24 mag; this 
results in an uncertainty of 0.039(1 + z). As a c r oss-ch eck, we 
also used the photometric redshifts of lllbert et al.1 d2009l) . which 
became available in the meantime. These are based on a photo- 
metric catalog of the COSMOS survey with 30 photometric bands 
including also narrow- and intermediate-band filters, and the ef- 
fects of both dust extinction and emission lines on template fitting 
were accounted for carefully during their computation. Based on 
the comparison with the spectroscopic redshift s of the / < 22.5 
mag zCOSMOS sample, the dispersion of the lllbert et~al1 (2009) 
photo-z catal og is ~0.007(1 + z) (and ~0.012(1 + z) down to 
/ = 24 mag; lllbert et alj|2009l) . Owin g to the use of a m uch larger 
number of photometric passbands, the lllbert et al.1 d2009h redshifts 
have a higher accuracy than our ZEBRA estimates. The advantage 
of the latter, however, is that they were derived by ourselves, and we 
are thus fully aware of their caveats and limitations. We therefore 
take the conservative strategy to use both versions of the COSMOS 
photo-zs in our analysis, as a means to identify any systematic un- 
certainty associated with the photometric redshift estimates. 

In order to estimate stellar masses from the photometric data, 
we used a newly-developed, non-public extension of our ZEBRA 
code, ZEBRA+, which estimates galaxy physical parameters, such 
as stellar masses, by fitting synthetic stella r population model 
SEDs to the galaxy broadband photometry (cf. iFontana et al.l2004l ; 
IPozzetti et al.ll2007l;lMarchesini et al1l2009l) . As in ZEBRA's Max- 
imum Likelihood mode, the best fit SED and its normalisation, are 
computed by minimising the x 2 between the observed and tem- 
plate fluxes. Additionally, ZEBRA+ includes reddening of the in- 
put SEDs with a variety of dust laws, when required. Using the 



1 The code is publicly available under the following URL: 

www . exp-astro .phys -ethz . ch/ ZEBRA 
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Figure 1. Example (HST ACS /-band) postage stamp images of both early- 
type (top two rows) and late-type (bottom two rows), barred and unbarred 
galaxies from our sample of 916 COSMOS disks at 0.2 < z < 0.6, M > 
10 10 ' 5 Mq, and e < 0.3. The width of each of these postage stamps was 
adjusted within a range of 3.3 to 6.3 arcsec in angular size to best emphasise 
the morphological details of each system. 



zCOSMOS 10k sample lBolzonella et all d2009h have demonstrated 
a level of uncertainty of aiogM ~ 0.20 arising from the choices 
of SED template and dust extinction law for stellar masses derived 
via SED-fitting to the COSMOS photometric data. This is consis- 
tent with the scatter of <7i og m ~ 0.27 between our ZEBRA+ stel- 
la r mass estimates and those from the publicly-available catalog 
of lPannella et alj d2009h for galaxies well-matched in both position 
(centroid difference less than 2 x 10~ 5 arcsec) and photometric red- 
shift (|zzebra - zparmdial < 0.02(1 + zzebra))- Details for ZE- 
BRA+, including a further investigation into the releva nt uncertain- 
ties, a re given in Oesch et al. (2010, in prep.; see also lOesch et al.l 
|2010|) . An acceptable synthetic template match for computation 
of stellar masses was able to be identified for all but 7033 of the 
109,446 galaxies with estimated photometric redshifts (a 6% fail- 
ure rate). We stress that our results remain unchanged when using 
either of the two photo-z catalogs; furthermore, the resulting un- 
certainties in the stellar mass and star formation rate estimates are 
negligible. A total of 32,840 galaxies were identified with photo- 
metric redshifts in the range 0.2 < z < 0.6 examined in this study, 
of which 3519 are above the adopted mass limit of M > 10 10 ' 5 

Mq. 

Star formation rates were computed based upon each galaxy's 
total IR (Lib.) luminosity (derived from the MIPS detections at 24 
um. lLeFloc'h et a l. 2009, d own to a flux limit of 0.15 mJy, using 
the templates of lRieke et al . 2009) and total UV (Luv) luminosity 
(estimated by interpolating the observed SEDs to 2800 A). Specifi- 



cally, we adopt the relation SFR = 9.8 x 10 _11 (-Lir + 2.2Ltjv) 
Mq yr" 1 from lBell etal] ( l2005h . The total uncertainties in our SFR 
estimates are ~0.4 dex, give n the ~0.3 dex intrinsic scatter about 
this relation dBell etall 2005) and the ~0.3 dex uncertainty on the 
estimation of UV luminosities from our SED template fits. Envi- 
ronmental densities are sourced from Scoville et al.'s (2007) cata- 



log of (large-scale) surface densities for COSMOS galaxies in fixed 
redshift intervals. 

Galaxy morphologies were obtained with an up graded version 
of the Zurich Estimator of Structural Types (ZEST; [Scarlata et al.l 
|2007|) . known as ZEST+ (Carollo et al. 2010, in prep). Relative to 
its predecessor, ZEST+ includes additional measurements of non- 
parametric morphological indices for characterising both structure 
and substructure. Moreover, ZEST+ offers a Support Vector Ma- 
chine (SVM) classification (adopted herein), as well as a Principal 
Component Analysis (PCA) scheme for consistency with the ear- 
lier ZEST. The ZEST+ SVM classifies galaxies in seven morpho- 
logical types located in specific regions of the 6-dimensional space 
of concentration, asymmetry, dumpiness, M20, Gini coefficient, 
and ellipticity. The different types were visually inspected to ensure 
a broad equivalence with the following Hubble types: E, SO-Sa, Sb, 
Sc, Sd, and Irr (with the Irr types divided into two classes, 'concen- 
trated' and 'non-concentrated', according to the value of the con- 
centration index). As a result of visual inspection of the 35 19 galax- 
ies in our sample at 0.2 < z < 0.6 and M > 10 10 ' 5 M Q , a total 
of 36 Irr, 14 SO-Sa, and 33 spiral (Sa-Sb to Sc-Sd) systems were 
identified as 'catastrophic' misclassifications and were manually 
reassigned to their respective types. The two most frequent causes 
of such 'catastropic' misclassifications were the presence of bright, 
overl apping neighbours and/or strong dust lanes (see lSargent et al.l 
120101 for a detailed study of the impact of dust on disks in COS- 
MOS), which can bias recovery of the quantitative morphological 
indices used as input to the SVM analysis module. Regarding the 
reliability of ZEST+ output in general, calibrations performed by 
Scarla ta et al. 1 mi) using the original ZEST code to classify real 
COSMOS galaxies show that down to I = 24 mag the misclassifi- 
cation rate is at most ~30% for ellipticals incorrectly identified as 
early- and intermediate-type disks in the lowest S/N bin. As ZEST+ 
features substantially improved algorithms for computation of the 
relevant quantitative morphological indices relative to the original 
ZEST, we are confident that the ZEST+ morphological classifica- 
tions are even more reliable than those of ZEST. 

The end result of our classification process was a final sample 
of 2820 disk galaxies (morphological type between SO-Sa and Sc- 
Sd), of which 916 have ellipticities (e) less than 0.3 (equivalent 
to disk inclination less than 45 deg). Example (HST ACS J-band) 
postage stamp images of both early-type and late-type, barred and 
unbarred systems from our final sample of 916 COSMOS disks at 
at 0.2 < z < 0.6, M > 10 10 ' 5 M Q , and e < 0.3 are presented in 
Fig. [2 



3 METHOD 
3.1 Bar Detection 

Barred galaxies in our sample were identified as follows. First, 
the ellipse package in I RAF was used to recover the ellipti- 
cal isophote profile of each galaxy from its ACS J-band image. 
Specifically, ellipse was run repeatedly for each galaxy, vary- 
ing the initial guesses of object position, semi-major axis, elliptic- 
ity, and position angle until a robust fit was recovered down to the 
level of sky noise (cf. Regan & Elm egreen 1 997l ; lLaine et al.l2002l ; 



iMenendez-Delmestre et al.ll2007l) . If a robust fit could not be iden 
tified with a single combination of initial parameter guesses (<5% 
of cases), the profile was constructed piece-wise from the robust 
sections of two or more fitting attempts. Pixels within 2.5 times 
the Kron radius of any nearby neighbours were excluded from the 
fitting process. 
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All galaxies with profiles displaying a monotonic increase 
in ellipticity of e m - m ^ 0.4 with a subsequent drop of at least 
Ae ^0.1 were flagged as candidate barred systems. The threshold 
ofe m i n ^ 0. 4 was chosen to rep licat e the selection of str ong galaxy 
bars made bv ljogee et al. | J2004I) and lSheth et al] J2008h . However, 
unlike these other studies, we did not reject candidate barred sys- 
tems based on a further criterion of isophote position angle change 
at the bar end. Rather, the final selection of barred galaxies was per- 
formed via visual inspection of the candidates, thereby maximising 
the completeness of our bar sample (i.e., allowing for detection of 
those galaxies in which, by chance, the uncorrelated bar and disk 
position angles are closely aligned; cf. Menendez-Delmestre et al. 
2007). In total, 213 (67%) of our 320 candidate barred systems 
(23% of the entire disk galaxy sample) were thereby confirmed as 
barred. A further 23 barred candidates that were ultimately rejected 
were noted as 'ambiguous' cases (e.g. the 'bar' morphology was 
perhaps more consistent with that of a lens or truncated inner disk). 
Example (HST ACS /-band) postage stamp images of both early- 
type and late-type, barred and unbarred galaxies from our final sam- 
ple of 916 COSMOS disks at at 0.2 < z < 0.6, M > 10 10 5 M©, 
and e < 0.3 are presented in Fig.Q] 

3.2 Accounting for Selection Effects 

Completeness limits on bar detection given the isophotal elliptic- 
ity criteria described above were estimated via artificial galaxy 
simulations. Both early-type and late-type barred galaxy models 
were constructed from three components: a Sersic bulge, a (trun- 
cated) Sersic bar, and an exponential disk. The properties of each 
model were chosen to be consistent with t he observed properties 
of barred galaxies in the local universe telmegreen et al.l l 19961 ; 
lLaurikainen et alj|2007l ; iGadottil |2008| ; IWeinzirl et al.l 120091) 7 and 
are listed in Table Q] The use of three different early type models 
was necessary to replicate the variety of structural types o bserved 
amongst local SO and Sa systems (e.g. iBalcells et al.l2003l) . Impor- 
tantly, the models adopted span a wide range of bulge Sersic indices 
(1 < n < 4) and bulge-to-total flux ratios (0.1 B/T 0.6), and 
thus encompass the entire parameter space observed for massive 
local disk galaxies. 

For all models we adopt a bar ellipticity of e = 0.6 as this 
represents a lower limit on the intrinsic bar ellipticities typically 
recovered for massive disk galaxies in 2D bulge-bar-disk structural 
decomposition studies (e.g. Gadotti 2008, MNRAS). The lower the 
intrinsic ellipticity of the bar, the lower the mean surface brightness 
(at a given Bar/T), so it follows that that low ellipticity bars will 
be most sensitive to surface brightness-dependent detection biases. 
Moreover, the lower the intrinsic bar ellipticity, the more sensitive 
will the system be to size-dependent detection biases caused by the 
'rounding' of bar isophotes by the PSF. Hence, by using a 'lower 
limit' bar ellipticity of e = 0.6 in all models we best expose the 
limits of our detection routine, and thereby identify safely conser- 
vative selection limits. 

Artificial galaxy images for each type of model were gener- 
ated and convolved with a Gaussian PSF of 0.12 arcsec (4 pixels) 
FWHM, then degraded to the noise level of the COSMOS ACS I- 
band imaging. A total of 1960 images were constructed for each 
model in a grid of total galaxy magnitudes and disk scalelengths 
spanning I = 17 to 24 mag and h = 0.15 to 3 arcsec (5 to 
100 pixels), respectively. Ten instances of each model were gen- 
erated at each grid point with random disk inclinations ranging to 
45 deg, and random component position angles (with the bar size 
scaled relative to the projection of the disk). These images were 



Table 1. Key Model Parameters for Artificial Galaxy Simulations of Bar 
Detection Completeness 



# 


Type 


B/T 


Bar/T 


Re/h a 


n a 


^bar 


L ba . T /h a 


I 


Early 


0.4 


0.15 


0.2 


4.0 


0.7 


1.15 


2 


Early 


0.2 


0.15 


0.2 


2.5 


0.7 


1.15 


3 


Early 


0.6 


0.15 


0.2 


4.0 


0.7 


1.15 


4 


Late 


0.1 


0.10 


0.1 


1.0 


0.9 


0.90 



a These symbols refer to the effective radius (R e ), Sersic index (n), and 
truncation length (L) of the generalised 2D Sersic profile (cf. Aguerri 
et al. 2009; Graham & Driver 2005) employed for the model bulge and 
bar components, while h is the scalelength of the underlying exponential 
disk. In all cases we adopt ellipticities (e) of e = and e baI = 0.6, and 
shape parameters (c) of c = 2 (elliptical isophotes) and c bar = 2.4 (boxy 
isophotes), for the bulge and bar respectively. Note that in each artificial 
galaxy image generated the bar length is scaled relative to the projection of 
the disk. 

run through our candidate bar detection pipeline, and the results 
are plotted as a function of total galaxy magnitude and total half 
light radius (more easily comparable against observations than disk 
scalelength) in Fig. [2] 

These simulations indicate a high bar detection completeness 
for most galaxies in our sample, supposing a mix of early-type 
and late-type disks equivalent in structure to our models #1, #2, 
and #4 (see TableQJ. Although simulations with early-type model 
#3 revealed that our detection procedure is unlikely to discover 
bars in any severely bulge-dominated systems (B/T > 0.6), high- 
resolution imaging studies indicate a neglig ible large-scale bar frac- 
tion for such galax ies in any case (e.g. ILaurikaine n et al.l 1200*71 ; 
IWeinzirl et alj|2009l) . Comparing the distribution of ZEST+ classi- 
fied disks with M > 10 10 ' 5 M© and i < 45 deg at 0.2 < z < 0.6 
in the COSMOS field against the model completeness limits at low 
surface brightness (for models #1, #2, and #4, see Fig.|2j reveals 
that bar detection in our sample is not surface brightness limited. 
In particular, there are very few galaxies with observed magnitudes 
and half light radii approaching these limits, even after account- 
ing for the systematic biases at low signal-to-n oise inherent in the 
Kron -style measurements employed here (cf. ICameron & Driver] 
2009). However, the limit for bar detection in galaxies with small 
apparent half light radii (i? e ~ 0.33 arcsec) is significant and must 
be considered in our analysis. 

To ensure a uniform detection completeness of barred disks 
at all redshifts we exclude from our subsequent analysis all galax- 
ies smaller than the physical size represented by this apparent size 
limit at the highest redshifts explored here (i.e., R e ,nm = 2.2 kpc 
at z = 0.6). This cut removes 283 (~31%) of our inital sample of 
916 COSMOS disks, including 30 bars confidently detected in (pri- 
marily early-type) galaxies with R e ~ 1.5—2.2 kpc at z < 0.45. 
Hence, we can confirm that bars are observed to exist in compact 
disks at these intermediate redshifts, but note that the evolution of 
this population cannot be robustly constrained with the available 
data. Therefore, we emphasise that the results presented in this 
study concern intermediate-to-large disks (R e > 2.2 kpc) only. 



4 RESULTS 

In Fig. [3] we present the (strong) bar fraction of ZEST+ classified 
disk galaxies with M > 10 10 ' 5 M©, i < 45 deg, and R e > 2.2 
kpc in the COSMOS field in redshift bins of 0.2 < z < 0.3, 
0.3 < z < 0.4, 0.4 < z < 0.5, and 0.5 < z < 0.6. Two key 
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Figure 2. The recoverability of bars in the COSMOS ACS /-band imaging for early-type (left) and late-type (right) model disks at low inclination (i < 45 deg) 
as derived via artificial galaxy simulations. The detection completeness of input bars identified via our isophotal ellipse-fitting procedure in each bin of intrinsic 
(i.e., input) apparent magnitude and half light radius is indicated by the large grey-shaded circles on a scale from white representing 100% completeness to 
black representing 0% completeness. The effective limits on (intrinsic) apparent size and surface brightness for high detection completeness are highlighted 
by dashed, red lines. Blue arrows indicate the mean errors on the measurement of (Kron-style) galaxy magnitudes and sizes towards low surface brightnesses, 
and red, solid lines mark boundaries within which galaxies are no longer scattered in (due to measurement errors) from regions of low detection completeness. 
Finally, the orange datapoints reveal the observed distribution of ZEST+ classified disk galaxies with M > 10 10 5 Mq and i < 45 deg at 0.2 < z < 0.6 in 
the COSMOS field. Note that only the results for simulations corresponding to our models #1 and #4 (see Table[T} are displayed here; the results for model #2 
were nearly identical to those for model #1, while almost no bars were recovered for model #3 (see the discussion in Section l3~2t . 




t - ^ 1 1 1 1 1 r^~i 1 1 1 1 r 



0.0 0.1 0.2 0.3 0.4 0.5 0.60.0 0.1 0.2 0.3 0.4 0.5 0.60.0 0.1 0.2 0.3 0.4 0.5 0.6 
Redshift Redshift Redshift 

Figure 3. The redshift evolution of the bar fractions of ZEST+ classified disk galaxies with M > 10 10 ' 5 Mq, i < 45 deg, and R e > 2.2 kpc at 
0.2 < z < 0.6 in the COSMOS field subdivided by stellar mass, morphological type, SSFR, and environmental density. Galaxies at intermediate stellar 
mass (10 10 5 < M < fO 11 Mq) are shown on the top row, while those at high stellar mass (M > fO 11 Mq) are shown on the bottom row. The left 
column displays the bar fraction evolution by ZEST+ morphological type for early-type and late-type disks. The middle and right columns display the bar 
fraction evolution by SSFR and (large-scale) environmental density respectively for the lowest and highest thirds of the population in each mass and redshift 
bin. In each panel the thick red and blue lines indicate the bar fractions of the relevant populations, and the corresponding uncertainties are revealed by the 
cross-hatched areas. The black squares and error bars overlaid on the plots by density indicate the total disk galaxy bar fraction in each redshift bin. Note the 
lack of evolution at M > IO^Mq, and the net evolution below this mass scale. The local Universe bar fractions marked at z = 0.03 indicate the values 
measured in our low redshift (0.02 < z < 0.04) SDSS comparison sample selected and analysed as described in the Appendix. 
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Table 2. Bar Fractions by Mass, Morphology, SSFR, and Environmental Density in COSMOS at 0.2 < z < 0.6 



Redshift Stellar Mass™ 

0.2 < z < 0.3 Intermediate Early-Type 
0.2 < z < 0.3 Intermediate Late-Type 
0.2 < z < 0.3 High Early-Type 
0.2 < z < 0.3 High Late-Type 


Morphological Type Specific 

/ bar = 0.44 ± 0.08 Low SSFR 
/ bar = 0.25 ± 0.13 High SSFR 
/bar = 0.21 ± 0.11 Low SSFR 
/bar = 0.50 ± 0.35 High SSFR 


Star Formation Rate™ En 

/bar = 0.50 ± 0.11 High density 
/bar = 0.27 ± 0.11 Low density 
/bar = 0.20 ± 0.18 High density 
/bar = 0.33 ± 0.19 Low density 


vironmental Density™ 

/bar =0.44 ±0.11 
/bar =0.38 ±0.11 
/bar =0.16 ±0.15 
/bar = 0.40 ± 0.21 


0.3 < z < 0.4 Intermediate Early-Type 
0.3 < z < 0.4 Intermediate Late-Type 
0.3 < z < 0.4 High Early-Type 
0.3 < z < 0.4 High Late-Type 


/bar = 0.44 ± 0.05 Low SSFR 
/bar = 0.29 ± 0.06 High SSFR 
/bar = 0.23 ± 0.06 Low SSFR 
/bar = 0.57 ± 0.18 High SSFR 


/bar = 0.54 ± 0.07 High density 
/bar = 0.24 ± 0.06 Low density 
/bar = 0.11 ± 0.07 High density 
/bar = 0.44 ± 0.12 Low density 


/bar = 0.39 ± 0.07 
/bar = 0.40 ± 0.07 
/bar =0.28 ±0.10 
/bar =0.17 ±0.09 


0.4 < z < 0.5 Intermediate Early-Type 
0.4 < z < 0.5 Intermediate Late-Type 
0.4 < z < 0.5 High Early-Type 
0.4 < z < 0.5 High Late-Type 


/bar = 0.35 ± 0.07 Low SSFR 
/ bar = 0.14 ± 0.05 High SSFR 
/bar = 0.25 ± 0.07 Low SSFR 
/ bar = 0.40 ±0.15 High SSFR 


/bar = 0.31 ± 0.07 High density 
/bar = 0.08 ± 0.05 Low density 
/bar = 0.27 ± 0.11 High density 
/bar = 0.44 ± 0.12 Low density 


/bar = 0.29 ± 0.07 
/bar =0.17 ±0.06 
/bar =0.13 ±0.08 
/bar =0.46 ±0.13 


0.5 < z < 0.6 Intermediate Early-Type 
0.5 < z < 0.6 Intermediate Late-Type 
0.5 < z < 0.6 High Early-Type 
0.5 < z < 0.6 High Late-Type 


/bar = 0.38 ± 0.06 Low SSFR 
/bar = 0.11 ± 0.04 High SSFR 
/bar = 0.26 ± 0.06 Low SSFR 
/bar = 0.43 ± 0.19 High SSFR 


/bar = 0.34 ± 0.07 High density : / ba r = 0.16 ± 0.05 
/bar = 0.03 ± 0.03 Low density : / ba r = 0.34 ± 0.07 
/bar = 0.16 ± 0.08 High density : / bar = 0.37 ± 0.11 
/bar = 0.58 ± 0.11 Low density : / bar = 0.21 ± 0.09 



a Here we define intermediate stellar mass as 10 10 5 < M < 10 11 Mq and high stellar mass as M > 10 11 Mq, and we identify the low and high SSFR 
(or environmental density) populations as the lower and upper thirds of the SSFR (or environmental density) distribution in each mass and redshift bin (as 
explained in Section|4). The COSMOS disk sample examined is also restricted to systems with e < 0.3 and R e > 2.2 kpc as described in Section [T2l 



stellar mass regimes are explored: intermediate mass 

(10 10.5 < 

M < 10 11 M ), and high mass (M > 10 11 M s ). Within each 
mass regime we further explore the dependence of bar fraction on 
(ZEST+ classified) morphological type (early-type vs. late-type), 
specific star formation rate, and (large-scale) environmental den- 
sity. For the purposes of exploring the latter two redshift depen- 
dencies we identify the lower and upper thirds of the population in 
each parameter in each mass and redshift bin. We also note that the 
error bars shown in Fig. [3] reflect the Poisson uncertainties only; 
cosmic variance potentially contributes an additional ~13% uncer- 
tainty fo r populations of this nu mber density in the volume sam- 
pled (cf. lTrenti & Stiavellill2008h . The bar fractions and statistical 
uncertainties for each redshift, mass, morphology, SSFR, and en- 
vironmental density bin measured here (and shown in Fig. [3) are 
compiled in Table[2]for reference. 

At intermediate stellar mass the mean bar fraction of early- 
type disks (/bar = 0.42 ± 0.05) is much higher than that of late- 
type disks (/bar = 0.19 ± 0.03) over the redshift interval ex- 
plored. However, at high stellar mass this trend is reversed with 
early-type disks less frequently barred (/bar = 0.20 ± 0.04) 
than late-type disks (/bar = 0.43 ± 0.07). The bar fractions 
of galaxies with the lowest and highest specific star formation 
rates at both mass regimes closely follow the trends evident in 
the morphologically-selected early-type and late-type populations 
respectively — reflecting a close relationship between disk morphol- 
ogy and star formation rate at these epochs. The bar fractions by 
large-scale environmental density reveal no relationship between 
the frequency of bars and environment at fixe d mass. However, we 
caution that the environmental densities from lScoville et all J2007I) 
used here are based on photometric redshifts and averaged over 
>100 Mpc scales. Hence, we are unable to confirm a null depen- 
dence of the bar fraction on density down to group scales with this 
dataset. Finally, the total bar fraction is observed to build up by a 
factor of ~2 between z ~ 0.6 and z ~ 0.2 at intermediate stellar 
mass, but at high stellar mass a null evolution scenario is favoured. 

The local Universe bar fractions measured in a complemen- 



tary sample of SDSS disks at 0.02 < z < 0.04 described in the 
Appendix are also shown in Fig.[3]for comparison. At intermediate 
stellar mass these low redshift datapoints confirm that the bar frac- 
tion in early-type disks remains constant, or even exhibits a slight 
decline, from z ~ 0.2 to 2 ~ while assembly of the late- type 
barred disk population continues to the present day. At high stel- 
lar masses these local Universe results are consistent with the null 
evolution of the early-type bar fraction since z ~ 0.6 evident in 
the COSMOS dataset. At face value, the bar fraction in late-type 
disks at high mass in the SDSS at 0.02 < 2 < 0.04 is much lower 
than that at 0.2 < z < 0.6 measured in COSMOS; however, the 
statistical uncertainites are rather large due to the low number den- 
sity of this galaxy population. Indeed, there is a clear dependence 
of the total bar fraction on mass in both our low and high redshift 
samples, supporting our main conclusions. 



5 DISCUSSION 

Several observational studies indicate that the large/massive disk 
galaxy population has reached a mature state by redshift unity. 
These studies include the null ev olution in the disk galaxy mass- 
size relation (Barden et al. 2005) , the size function of large disk s 
dLillv et alj 1 19981 : ICameron & Driver! 12001 ISargent et all l2007h. 
the 5 -band Tully-Fisher relation for massive disks Jziegler et al~ 
2002), and the morphological mix of massive galaxies dOesch et~ 
20101) since z ~ 0.7—1, as well as the relati ve paucity of highly dis - 
turbed and highly diffuse galaxies at z < 1 dConselice et al. 2004). 
Our work adds an important additional piece of information: There 
exist substantial populations of barred early- and late-type disks — 
whether defined by morphology or SSFR — at M > 10 11 Mq out 
to 2 ~ 0.6. At such high stellar masses the fractions of barred 
disks in each of these galaxy classes remain roughly constant over 
the ~3.3 Gyr that separate the 2 ~ 0.6 from the 2 ~ 0.2 epochs. 
Thus, already at the intermediate redshifts explored in our study, 
spectral and morphological transformations within the most mas- 
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sive disk galaxies have largely converged, to such a level of detail, 
to the familiar Hubble sequence that we observe in the local Uni- 
verse. 

In contrast, at intermediate stellar mass (10 10,5 < M < 10 11 
Mq) there is a substantial build-up in the bar fractions of early- and 
late-type disks, and of the total disk population (by a factor ~2), 
over the past ~3.3 Gyr of cosmic time. At such intermediate mass 
scales, therefore, the disk galaxy population — and the disk-galaxy 
sector of the Hubble sequence — continues structural evolution until 
at least z ~ 0.2. The preference for (large-scale) bars in intermedi- 
ate mass_g_ahKiesJoreside in early-type disks — a trend also noted 
bv lSheth etall d2008l) in their analysis of the COSMOS field— is 
perhaps a reflection of their earlier formation times compared to 
late-type disks of similar mass, i.e., these galaxies have evolved for 
longer in a dynamical state conducive to bar formation. 

Interestingly, as the bar fraction reveals no dependence on 
large-scale environmental density at fixed mass, it would appear 
that any enhancement in th e minor merger rate with environment 
(e.g. lHeiderman et alj20 09) at these mass scales and redshi fts does 
not lead to significant disk heating (cf. Hopkins et al. 2008), which 
could strongly impede bar formation l lAth anassoula & SellwooJ 
1986). Although, an increase in the rate of tidally-triggered bar for- 
mation in dense environments could also be masking this effect. 
We stress, however, that the 'environment' in our analysis, based 
on photometric redshifts, is averaged over >100 Mpc scales. The 
study of the bar fraction as a function of the group environment at 
early epochs is an important open qu estion that we are ex ploring 
using the zCOSMOS group catalog of iKnobel et al.l d2009h . More- 
over, given the importance of the grou p environment in dete rmin- 
ing the properties of local galaxies (e.g. lRobotham et alj20ld) . we 
will also soon be using the ZENS dataset (Carollo et al., in prep.) 
to construct a key benchmark for the barred galaxy fraction in local 
groups (Cameron et al., in prep.). 

Finally, our study also reveals the importance of investigating 
the evolution of barred galaxies at fixed stellar mass: this provides 
evidence for a substantial change in the relative fractions of bars in 
early- and late-type disks, and in the evolution of the total bar frac- 
tion, when stellar mass straddles across the 'critical' mass scale 
of ~1O 11 M0. The close relationship between galaxy evolutionary 
histories and stellar mass has also been demonstrated recently via 
studie s of local ga laxy morphologies dBundv. Ellis. & Conselicej 
2005) an d colours (Baldrv et al. 2006) , as well as in h igh redshift 
samp les dKovac et alj l20ld : iBolzonella et al.l 1 20091 : iTasca et ail 
2009). We therefore highlight the importance of constructing even 
larger, mass-limited disk galaxy samples at low and high redshift 
for improving our understanding of bar formation — an endeav- 
our that will be much assisted by various ground-based (UKIDSS, 
VIKING) and space-based (GOODS-ERS, HUDF09) near-infrared 
imaging surveys now under way, and their successors. 
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APPENDIX A: LOW REDSHIFT COMPARISON SAMPLE 

To provide a local Universe benchmark for comparison against 
the high redshift bar fractions quantified here in the COSMOS 
dataset we have constructed a complementary sample of low red- 
shift disks drawn from the Sloan Digital Sky Survey (SDSSfl 
From the visually-cla ssified SDSS galaxy catalog presented by 
iNair & Abraham! j2010h we have selected a subsample of 65 1 disk 
galaxies (SO- to Sd, -3 < visual class < 7) at 0.02 < z < 0.04 
with M > 10 10,5 M , axial ratio > 0.7 (i.e., disk inclination 
< 45 deg), and R e > 2.2 kpc. At these redshifts the physical 
resolution scale (0.6-1.1 kpc/FWHM) and rest-frame wavelength 
coverage (~4000-7000 A) of the ground-based SDSS g and r-band 
imaging are roughly comparable to those of galaxies in our primary 
sample at 0.2 < z < 0.6 in the space-based COSMOS /-band 
imaging (~4300-7000 A and 0.4-0.8 kpc/FWHM). 

To ensure consistency with our selection of barred galaxies in 
COSMOS we have applied an identical ellipse-fitting bar detection 
procedure to the SDSS comparison sample. For this purpose we 
downloaded the u, g, r, i, and z -band images of all 65 1 disk galax- 
ies in our local Universe sample from the SDSS imaging server 
(http ://das.sdss.org/ imaging/). Following the process 
described in Section [3TTI we recovered elliptical isophote profiles 
for each galaxy in the g and r-bands using ellipse, and identi- 
fied (strong) barred galaxy candidates as those objects with maxi- 
mum ellipticity e 0.4 and a subsequent drop of Ae ^0.1. Visual 
inspection of the u, g, r, i, and z-band images of the barred galaxy 
candidates was then used to construct master catalogs of confident 
g and r-band bar detections. During the visual inspection process 6 
galaxies with axial ratio > 0.7 were noted to be bulge-dominated 
SOs with near edge-on (rather than face-on) disks and were thus 
excluded in the subsequent analysis. 

From the 588 intermediate mass (10 10 5 < M < 10 11 M a ) 
disk galaxies in our local Universe sample we identify 214 barred 
systems (a bar fraction of 36%) in the g-band and 240 barred sys- 
tems (41%) in the r-band. From the 57 high mass (M > 10 11 
Mq) disk galaxies we identify only 9 barred systems (16%) in the 
g-band and 10 barred systems (18%) in the r-band. The difference 
of ~12% between the bar fractions recovered here in the g and r- 
bands is slightly higher than that of ~5% reported in the analysis 
of a hi gh resolution {D < 100 Mpc) SDSS sample bv lSheth et ail 
(2008). However, we note that the majority of bars missed by our 
g-band ellipse-fitting lie in compact galaxies with 2.2 < R e < 3 
kpc. Hence, we suspect that the larger PSF FWHM in the SDSS g- 
band images (FWHM ~ 1.5 arcsec) relative to the r-band (FWHM 
~ 1.4 arcsec) was the cause of the bar detection failures for these 
compact systems. At R e > 3 kpc the difference in recovered bar 
fractions between the g and r-bands in our sample is only 6%. In 
Fig. lAll we contrast the dependence of bar fraction on disk morpho- 
logical type (at R e > 2.2 kpc) recovered using our g and r-band 
ellipse-fit catalogs in order to demonstrate the minimal dependence 
of bar detection efficiency on (optical band) wavelength across all 
morphological types. 

In Table I A 1 1 we present the bar fractions in our local Uni- 
verse comparison sample subdivided by mass, morphological type, 



Bars in early- and late-type disks in COSMOS 9 



SDSS Visual Morph. Cat. 0.02 < z < 0.04 




Figure Al. The bar fraction by visual morphological class for massive 
(M > 10 10 ' 5 Mq) disks at axial ratio > 0.7 and R e > 2.2 kpc in the lo- 
cal Universe (0.02 < z < 0.04). Th e visual morphological classifications 
are from the lNair&Abrahaml <2010l) catalog. The bar fractions recovered 
via application of the ellipse-fitting method (described in Section l3Tt to the 
SDSS g- and r-band images of each galaxy are compared here to demon- 
strate the minimal dependence of bar detection efficiency on (optical band) 
wavelength across all morphological types. 

SSFR, and environmental density as performed for our high 
redshift disk sample in Section [4] Specifically, we employ the 
iKauffmann et al] J2003h stellar masses to separate our SDSS sam- 
ple into intermediate mass systems at 10 10 ' 5 < M < 10 11 Mq 
and high mass systems at M > 10 Mq. Early-type and late-type 
disks are defined as those with visual morphological classes SO- to 
Sa and Sb to Sd, respectively. To define our low (high) SS FR and 
density popul ations we adopt the iBrinchmann et al] l l2004r) SSFR 
estimates and iBaldrv et al] 120061) environmental densities (aver- 
age of S4 and E5), and identify the lowest (highest) thirds of each 
distribution. 

A comparison of these local Universe bar fractions against 
those recovered at high redshift in COSMOS is presented in Sec- 
tion|4] As described therein, the SDSS bar fractions by stellar mass 
and spectral/morphological type confirm the evolutionary trends 
observed for massive barred galaxies over the redshift interval 
0.2 < z < 0.6 in COSMOS. 



2 Whilst there already exists a valuable database of SDSS systems artifi- 
cially redshifted to 0.7 < 2 < 1.2 and matched to the P SF and signal-to- 
noise of the COSMOS ACS I-band imaging provided by Kampc zvk et all 
120071) . the effective resolution of these artificially redshifted galaxies is 
much lower than that for the 0.2 < z < 0.6 COSMOS disk sample em- 
ployed in this study, and thus not suitable for the present analysis. 
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Table Al. Bar Fractions by Mass, Morphology, SSFR, and Environmental Density in an SDSS Local Universe Comparison Sample 



Bar Catalog Stellar Mass" Morphological Type Specific Star Formation Rate" Environmental Density" 



r-band ellipse fit 


Intermediate 


Early-Type : 


/bar 


= 0.40 ± 0.03 


Low SSFR : 


/bar 


= 0.49 ± 0.04 


High density 


■ /bar 


= 0.48 ± 0.04 


r-band ellipse fit 


Intermediate 


Late-Type : 


■ /bar 


= 0.41 ±0.03 


High SSFR : 


■ /bar 


= 0.34 ±0.03 


Low density 


• /bar 


= 0.39 ±0.03 


r-band ellipse fit 


High 


Early-Type : 


■ /bar 


= 0.13 ±0.06 


Low SSFR : 


/bar 


= 0.21 ±0.09 


High density 


• /bar 


= 0.21 ±0.09 


r-band ellipse fit 


High 


Late-Type : 


■ /bar 


= 0.22 ±0.08 


High SSFR : 


/bar 


= 0.05 ± 0.05 


Low density 


• /bar 


= 0.16 ±0.08 



a Here we define intermediate stellar mass as 10 10,5 < M < 10 11 M and high stellar mass as M > 10 11 Mq, and we identify the low and high SSFR 
(or environmental density) populations as the lower and upper thirds of the SSFR (or environmental density) distribution in each mass and redshift bin. The 
SDSS disk sample examined is also restricted to systems with axial ratio > 0.7 and R e > 2.2 kpc as described in this Appendix. 



© 0000 RAS, MNRAS 000, 000-000 



